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Photoacoustic imaging is a novel hybrid imaging modaUty combining the high spatial resolution of optical 
imaging with the high penetration depth of ultrasound imaging. Here, for the first time, we evaluate the 
efficacy of various photosensitizers that are widely used as photodynamic therapeutic (PDT) agents as 
photoacoustic contrast agents. Photoacoustic imaging of photosensitizers exhibits advantages over 
fluorescence imaging, which is prone to photobleaching and autofluorescence interference. In this work, we 
examined the photoacoustic activity of 5 photosensitizers: zinc phthalocyanine, protoporphyrin IX, 2,4-bis 
[4-(N,N-dibenzylamino)-2,6-dihydroxyphenyl] squaraine, chlorin e6 and methylene blue in phantoms, 
among which zinc phthalocyanine showed the highest photoacoustic activity. Subsequently, we evaluated its 
tumor localization efficiency and biodistribution at multiple time points in a murine model using 
photoacoustic imaging. We observed that the probe localized at the tumor within 10 minutes post injection, 
reaching peak accumulation around 1 hour and was cleared within 24 hours, thus, demonstrating the 
potential of photosensitizers as photoacoustic imaging contrast agents in vivo. This means that the known 
advantages of photosensitizers such as preferential tumor uptake and PDT efficacy can be combined with 
photoacoustic imaging capabilities to achieve longitudinal monitoring of cancer progression and therapy in 
vivo. 



Photoacoustic (PA) imaging is a rapid emerging biomedical imaging modality which provides in vivo 
functional imaging information at clinically relevant penetration depths, while maintaining high spatial 
resolution and image contrast, as compared to existing imaging techniques' ^ In PA imaging, the tissue is 
irradiated with a laser pulse and the light energy is absorbed by endogenous chromophores such as hemoglobin 
and melanin, as well as exogenous contrast agents in tissue. This causes the tissue to undergo thermoelastic 
expansion and generate corresponding acoustic pressure waves, which in turn can be detected via ultrasound 
transducer arrays. The higher penetration depth of PA imaging (5-6 cm)'"^ over fluorescence and optical 
coherence tomography (OCT)^ enables deep tissue imaging, especially in clinical settings\ 

In terms of applications, endogenous hemoglobin in blood has been used for PA imaging of the tumor vascular 
network in the rat brain^, the blood- oxygenation dynamics in the mouse brain^'^, the human arm^, as well as breast 
imaging^. In addition, various exogenous contrast agents have been introduced to enhance the imaging con- 
trast'°"^\ such as carbon nanotubes (SWNTs)'^"'^, near-infrared (NIR) dyes like indocyanine green (ICG)'^"'^, as 
well as gold nanoparticles'^"^\ However, the clinical application of these contrast agents has been limited due to 
cytotoxicity issues. There have been various attempts to circumvent this problem, such as the surface modifica- 
tion of gold nanorods with polyethylene glycol (PEG)^^, in order to lower the cytotoxicity and increase the blood 
circulation time, but there still remains a need for a great leap forward towards clinical translation. 

On the other hand, photodynamic therapeutic (PDT) agents, also known as photosensitizers, have been widely 
used in clinical trials for fluorescence imaging and PDT^^. These photosensitizers are generally classified as 
porphyrins or nonporphyrins^^. Porphyrin-based photosensitizers are further categorized as first-, second- or 
third- generation photosensitizers, in which protoporphyrins, chlorins and phthalocyanines fall under the sec- 
ond-generation category, whereas squaraines and methylene blue form part of the nonporphyrin-based photo- 
sensitizers. In this work, we examined the potential of five photosensitizers as PA contrast agents namely zinc 
phthalocyanine (ZnPc), protoporphyrin IX (PpIX), 2,4-bis [4-(N,N-dibenzylamino)-2,6-dihydroxyphenyl] 



SCIENTIFIC REPORT; j 4 : 5342 j DO!: 1 0.1 038/srep05342 



1 




700 750 800 850 900 

Wavelength (nm) 



Figure 1 | Normalized absorbance as a function of wavelength for the five 
photosensitizers in the NIR region from 680 to 900 nm. ZnPc shows the 
highest overall absorbance, followed by MB, Sq, PpIX and Ce6, in 
decreasing order. 

squaraine (Sq), chlorin e6 (Ce6) and methylene blue (MB), which, 
though not exhaustive, are representative of the various categories of 
photosensitizers, many of which have been used in clinical trials^^.To 
the best of our knowledge, this is the first evaluation of the various 
classes of photosensitizers as potential PA contrast agents. These 
photosensitizers exhibit low fluorescence quantum yields and thus 
can potentially possess high PA activity, since an excited system can 
either relax back to the ground state through fluorescence or ther- 
mally through internal conversion^^. Moreover, PA imaging of 
photosensitizers exhibits advantages over fluorescence imaging, 
which is prone to photobleaching and autofluorescence interference. 
Because these PDT agents preferentially accumulate in tumor due to 
the enhanced permeation and retention (EPR) effect^^, they offer 
tumor-targeted PA imaging. In this paper, we first evaluated the 
PA performance of the five PDT agents in a scattering phantom. 
Next, in order to demonstrate in vivo PA activity and tumor-target- 
ing efficacy, we injected ZnPc intravenously into mice and moni- 
tored the biodistribution over time using PA imaging. In light of the 
unique advantages such as clinical relevance, passive tumor-target- 
ing ability and high PA activity, these photosensitizer-based PA con- 
trast agents offer great potential in cancer diagnosis and therapy. 

Results and discussion 

In order to obtain information about the optical absorption prop- 
erties of the contrast agents, their wavelength-dependent absorption 
spectra at known concentrations were measured using a spectropho- 
tometer, which in turn were normalized to obtain the concentration - 
independent normalized absorbance (directly proportional to 
extinction coefficient) as a function of wavelength for each contrast 
agent (Figure 1), based on Beer-Lambert law: 

A{X) = ln{j^=f^{X)cl 

where A(k)is the measured wavelength- dependent absorbance, Iq is 
the incident light intensity, 1 is the transmitted light intensity, ef/lj is 
the wavelength-dependent extinction coefficient of the contrast 
agent, c is the concentration of the contrast agent and / is the path 
length of the quartz cuvette (1 cm). 

As shown in Figure 1, all the studied contrast agents exhibit a 
general decreasing trend in absorbance with increasing wavelength 
in the near infrared (NIR) range from 680 to 900 nm, in which ZnPc 
shows the highest overall absorbance, followed by MB, Sq, PpIX and 
Ce6 in decreasing order. Although these probes have higher reported 
molar extinction coefficients in the visible light range of the electro- 



magnetic spectrum^^, and hence, are expected to provide stronger PA 
signals at these wavelengths, we have chosen to study the PA activity 
of these compounds in the NIR region from 680 to 900 nm, because 
at these wavelengths there is lower tissue absorption enabling 
imaging at greater depths. At wavelengths above 900 nm, water 
absorption increases significantly, making the 680-900 nm range 
optimal for in vivo deep tissue imaging. 

In order to verify the PA activity of these photosensitizers under 
controlled conditions, we performed phantom measurements for all 
the contrast agents, in the 680-900 nm wavelength range, each at 5 
different concentrations. The phantom is cylindrical with a diameter 
of 2 cm and contains 2 cylindrical channels in which contrast agents 
can be placed to measure the PA signal, compared to a control agent 
(see Figure 2A). During the data acquisition, we recorded data from 
multiple transverse slices across the channel portion which contains 
the probe and control, and applied an excitation wavelength scan 
from 680 to 900 nm with an interval of 10 nm for each transverse 
slice, and recorded the averaged PA signals from 10 frames for each 
wavelength and position. After image reconstruction, results showed 
that all the contrast agents exhibit wavelength-dependent PA activity 
in phantoms. In addition, there is a similar trend in waveform 
between absorbance and PA intensity (both normalized) as a func- 
tion of wavelength for all the contrast agents, with a peak at around 
680-700 nm for both absorbance and PA intensity, that tapers 
downwards towards longer wavelengths. This demonstrates a strong 
correlation between optical absorption spectra and the PA spectra, 
which vahdates our data, as shown in Figure 2(B-F). In PA imaging, 
the PA intensity induced by optical absorption is proportional to 
light energy deposition, which is the product of the absorption coef- 
ficient and the local light fluence. Thus, the small deviations in trend 
between absorbance and PA intensities can be attributed to light 
fluence changes caused by slight variations in laser intensity. 

In addition, the PA signal for each contrast agent was spectrally 
unmixed via linear regression. This allows the isolation of the indi- 
vidual contribution of the contrast agent of interest that can be 
plotted as a function of concentration, which in turn was used to 
produce a straight line of best fit based on least- squares regression, 
for each contrast agent. As shown in Figure 3, the line corresponding 
to ZnPc has the highest gradient, which corresponds to the highest 
increase in PA signal for an incremental increase in concentration, 
when compared to the rest. We hereby define this gradient as a form 
of relative PA quantum yield (cpP), a kind of measure of the efficiency 
of the conversion of light absorption into PA signal. This is analogous 
to the fluorescence quantum yield (cpF), which measures the effi- 
ciency of the conversion of light absorption into fluorescence emis- 
sion. The relative cpP of the 5 compounds are computed and listed in 
Table 1. These are not absolute values, but arbitrary ratios, which 
reflect the relative PA strength of one compound against that of 
another. As shown in Table 1, it has also been reported that these 
photosensitizers have low fluorescence quantum yields (<0.2) but 
reasonably high singlet oxygen quantum yields (—0.5). This means 
that although these photosensitizers can offer high PDT efficacy, they 
exhibit low fluorescence, which may not be adequate for high con- 
trast diagnostic imaging. However, this limitation can be circum- 
vented by using PA imaging, as shown in this work. 

Moreover, on testing these compounds for dark toxicity on oral 
squamous carcinoma cell line (OSCC) for varying concentrations (1 
to 100 jiM), they exhibited more than 50% cell viability for concen- 
trations up to 100 laM of PpIX and MB, 50 |iM of Ce6 and 10 |iM of 
Sq and ZnPc (supplementary information. Fig SI). This vast range of 
concentrations tested was inclusive of the in vivo dose range for PDT 
for most of them (see SI Table SI) thus making a point that they can 
be used as PA imaging contrast agents at concentrations that do not 
induce much toxicity in vivo. 

Based on the prioritization of the different photosensitizers for 
MSOT imaging by phantom analysis, we next evaluated ZnPc in vivo. 
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Figure 2 | (A) Phantom schematic. (B-F) Normalized absorbance (black line) and PA intensity (colored line) as a function of wavelength. The strong 
correlation between optical absorption and the PA effect is illustrated by the similarity in trend between the two graphs for each contrast agent, with 
a peak at around 680-700 nm, that tapers downwards towards longer wavelengths. 
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Figure 3 | Multispectrally unmixed PA signal intensity as a function of 
concentration with a line of best fit for each contrast agent. 



since it was found to have the strongest PA signal in phantoms. In 
essence, we injected ZnPc intravenously into mice bearing subcutan- 
eous tumor and monitored the probe uptake over time in the liver, 
spleen, kidneys, intestines and tumor using PA imaging. We scanned 
a volume ROI consisting of transverse slices spanning from the liver 
to the lower abdomen and constructed maximum intensity projec- 
tions (MIPs) based on these slices for image analysis. A strong in vivo 
MSOT signal of ZnPc was detected in the reticuloendothelial system 
(liver, spleen), intestines and tumor site within the first hour after 
injection, with peak accumulation at the 1-hour time point, as shown 
in the background-corrected images in Figure 4. In addition, the 
non-background-corrected images in Figure 5 and Figure 6A shows 
that this signal gradually drops in intensity beyond 1 hour up to a 
day, suggesting probe clearance within a day. This allows for the 
rapid evaluation of PDT agent delivery after administration and 
can aid in optimizing PDT and not full form time point planning. 
Apart from monitoring the accumulation in tumor, the whole-body 
biodistribution of compounds can be visualized and quantified 
simultaneously by MSOT, which can aid in the elucidation of clear- 
ance pathways and choosing of optimal dosing strategies. We hypo- 
thesize that ZnPc is most likely cleared through the hepatobiliary 
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Table 1 Fluorescence, singlet oxygen and relative PA quantum yields of various photosensitizers 
Photosensitizers Fluorescence Quantum Yield (cpF) Singlet Oxygen Quantum Yield (cpA) 


Relative PA Quantum yield (cpP) 


ZnPc 0.06 [23] 


0.62 [23] 


0.47 


PplX 0.16 [23] 


0.56 [23] 


0.093 


Sq 0.037 [27] 


0.61-0.74 [28] 


0.18 


Ce6 0.19 [23] 


0.65 [23] 


0.036 


MB 0.02 [23] 


0.55 [23] 


0.029 



system, because high concentrations are found in the Hver and intes- 
tines up to the 3 -hour time point, whereas renal signal is much 
weaker in comparison, as shown in Figure 5 and Figure 6A. In addi- 
tion. Figure 6B also shows a peak tumor-to-muscle ratio of the 
MSOT signal at the 1-hour time point, making this time point ideal 
for PDT illumination and subsequent monitoring via MSOT 
imaging. Thus, in summary, the strong MSOT signals from ZnPc 
observed in the tumor region at the 1-hour time point suggests that 
MSOT imaging, combined with PDT, offers a novel theranostic 
approach with high translational potential. 

Conclusion 

In this work, we demonstrated that various classes of photosensiti- 
zers can be used as efficient novel contrast agents for PA imaging. 
Initially, we carried out the PA study in a phantom model and com- 
pared the performance of various classes of photosensitizers. Among 
the tested compounds, ZnPc was found to be the strongest in terms of 
PA signal with the highest relative PA quantum yield. Subsequently, 
we studied the in vivo localization and biodistribution of ZnPc in a 
preclinical animal model of cancer. The preferential accumulation of 
photosensitizers in tumors is due to the enhanced permeation and 
retention (EPR) effect that enables us to achieve longitudinal mon- 
itoring of cancer. In addition, unlike fluorescence imaging, PA 
imaging of photosensitizers is free from interference from photo - 
bleaching and autofluorescence. Moreover, MSOT imaging allows 
us to detect endogenous absorbers such as oxygenated and deoxyge- 
nated hemoglobin so that tumor blood vasculature and oxygenation 
status can be monitored at the same time. In addition, it has been 
reported that many well-studied photosensitizers have low fluor- 
escence quantum yields (<0.2) but reasonably high singlet oxygen 
quantum yields (—0.5). This means that although these photosensi- 
tizers can offer high PDT efficacy, they exhibit low fluorescence, 
which may not be adequate for high contrast diagnostic imaging. 
In this study, we have successfully demonstrated that these 




Figure 4 | (A) In vivo background-corrected MIP images of transverse slices through mouse preinjection and within the first hour after injection, 
showing gradual probe accumulation within the tumor site and various organs over time. (B) Time-resolved color coding of the MSOT signal within the 
first hour after injection, demonstrating peak localization of probe at the tumor site and various organs at 1 hour postinjection. A stack of MIPs of ZnPc 
signal over time (T = 10-60 mins) is condensed in a parametric map visualizing the T^ax of ZnPc; peak concentrations of ZnPc are observed at 
T = 60 mins. 



compounds offer reasonably strong PA signals, which offers a poten- 
tial alternative imaging technique, on top of fluorescence imaging. In 
future studies, more work will be devoted towards in vivo multifunc- 
tional imaging and therapeutic initiatives, such as longitudinal PA 
monitoring of tumor regression in response to photodynamic ther- 
apy for instance, as well as the design of drug carrier systems such as 
liposomes^^"^^ for reduced probe aggregation in blood, as well as 
enhanced tumor delivery and uptake. 

Methods 

Preparation of PA contrast agents. Zinc phthalocyanine, protoporphyrin IX, 
methylene Blue and 2,4-bis [4-(N,N-dibenzylamino)-2,6-dihydroxyphenyl] 
squaraine in powder form were purchased from Sigma- Aldrich. Chlorin e6 was 
purchased from Frontier Scientific (Logan, UT). Stock solutions of all compounds 
were prepared in DMSO. Formulations for absorbance measurements and animal 
studies were made in 10% DMSO. 

UV Absorbance measurement. The wavelength -dependent absorption spectrum of 
the contrast agents was measured from 680 to 900 nm using a spectrophotometer 
(DU 730, Beckman Coulter), which in turn was used as an input spectrum to be used 
for multispectral post-processing of PA signals. 

MSOT experimental parameters and protocol. All phantom and in vivo mouse 
imaging experiments were performed using a real-time multispectral optacoustic 
tomographic (MSOT) imaging system; inVision 128 (iThera Medical GmbH, 
Neuherberg, Germany). The phantom is made of polyurethane, cylindrical in shape 
with a diameter of 2 cm, which is specially designed to mimic the shape, size and 
optical properties of the mouse (iThera Medical GmbH, Neuherberg, Germany). In 
addition, it has 2 inner cylindrical channels, each with a diameter of 3 mm, one for 
holding the control medium and the other for holding the dissolved contrast agent in 
the same medium, as shown in Figure 2A. Optical excitation was provided by an 
optical parametric oscillator (OPO) with a tunable NIR wavelength range from 680 to 
980 nm, which is in turn pumped by a Q-switched Nd:YAG laser with a pulse 
duration of 10 ns and repetition rate of 10 Hz. Light was delivered by a fiber bundle 
divided into 10 output arms to illuminate the sample from multiple angles around the 
imaging plane. PA signals were acquired using a 128-element concave transducer 
array spanning a circular arc of 270°. This transducer array has a central frequency of 
5 MHz, which provided a transverse spatial resolution in the range of 150-200 |im. 
One transverse image slice was acquired from each laser pulse, resulting in a 
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Figure 5 | In vivo non-background-corrected MIP images of transverse 
slices through mouse at various time points postinfection, 
demonstrating strong probe MSOT signals within the liver, spleen, 
kidneys, intestines and tumor site. There is a gradual decrease in the tumor 
MSOT signal over time, which in turn implies probe clearance within a 
day. 

frame-rate of 10 Hz. During image acquisition, the sample is translated through the 
transducer array along its axis across the volume ROl, in order to capture the 
corresponding transverse image slices. For in vivo imaging, ultrasound gel was 
applied on the mouse skin surface and measurements were recorded in a 
temperature-controlled water for good acoustic coupling. An animal holder with a 
thin polyethylene membrane was used to prevent direct contact between the mouse 
and the water. 

Animal preparation. All animal experimental procedures were performed in 
accordance with the protocol #120774 approved by the Institutional Animal Care and 



Use Committee (lACUC). Xenograft mice models were established by injecting 
subcutaneously into the right flank of mouse, 0.2 mL of cell suspension containing 5- 
6 X 10*^ MCF-7, a human metastatic breast cancer cell line or OSCC, a human oral 
squamous carcinoma cell line and matrigel (BD biosciences) in 1 : 1 volume ratio. 
When the tumor volume reached a palpable size, the mouse was used for in vivo PA 
imaging. 

In vivo longitudinal monitoring of probe biodistribution in mouse xenograft 
model. ZnPc at a dosage of 0.21 mg/kg was injected via a catheter into the tail vein of 
tumor-bearing mice anaesthetized under isoflurane and the probe biodistribution 
was monitored over time in various organs using PA imaging. Before image 
acquisition, a volume ROI consisting of transverse slices with a step size of 0.3 mm 
spanning from the liver to the lower abdomen was selected by manual inspection of 
live MSOT images, and the 6 laser excitation wavelengths of 680, 700, 750, 800, 850 
and 900 nm were selected for correspondence with the major turning points in the 
absorption spectra of ZnPc, oxy- and deoxy-haemoglobin. Multispectral imaging was 
then performed with 10 signal averages per wavelength per transverse slice, before 
injection, during injection and 1, 3, 5, and 24 hr postinjection. 

Image reconstruction and multispectral processing. Images were reconstructed 
using a model-based approach^^ for offline analysis. After image reconstruction, 
spectral unmixing was performed to resolve individual components from different 
chromophores in the system. For each pixel in the image, the method fits the total 
measured optoacoustic spectrum to the known absorption spectra of the individual 
chromophores, based on least- squares linear regression. 

Image processing. Maximum intensity projection (MIP) images were prepared and 
presented in Figure 4 and 5, for better display of anatomy and quantification. For the 
first hour in Figure 4, the difference between the time-point images and the prescan 
image (before injection) were displayed, which is possible because the animal is intact 
and anaesthetized in the MSOT machine within the 1st hour. For the later time points, 
this background subtraction is not possible, as the animal is removed and placed in 
the MSOT system repeatedly, with repositioning, in order to avoid overdosage of 
isoflurane for the entire 24-hour duration. Although Figure 5 is not background- 
corrected, it still shows the biodistribution trend over the entire duration. 




10 20 30 40 50 60 180 300 1440 
Time post injection (min) 

Figure 6 | (A) Biodistribution of ZnPc at various organs and tumor at points. (B) Tumor-to-muscle ratios of MSOT signal at multiple time points (n = 
4). Tumor-to-muscle ratio of the MSOT signal peaks at 1-hour, thus making this time point ideal for PDT illumination and subsequent monitoring via 
MSOT imaging. 
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